Lattice defects and boundaries in conducting carbon nanotubes 
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We consider the effect of various defects and boundary structures on the low energy electronic 
properties in conducting zigzag and armchair carbon nanotubes. The tight binding model of the 
conduction bands is mapped exactly onto simple lattice models consisting of two uncoupled parallel 
chains. Imperfections such as impurities, structural defects or caps can be easily included into 
the effective lattice models, allowing a detailed physical interpretation of their consequences. The 
method is quite general and can be used to study a wide range of possible imperfections in carbon 
nanotubes. We obtain the electron density patterns expected from a scanning tunneling microscopy 
experiment for half fullerene caps and typical impurities in the bulk of a tube, namely the Stone- 
Wales defect and a single vacancy. 



I. INTRODUCTION 

The quasi one dimensional nature and interesting 
electronic properties of single-walled carbon nanotubes 
(SWCNTs) have attracted a lot of attention ever since 
they where first synthesized in 1993jii£ The tubes are well 
described by a graphene sheet that is rolled up along a 
wrapping vector (C/^ = A^ai + Ma2) with characteristic 
chiral indices (A/", M).^ The finite circumference results in 
discrete one-dimensional bands, which may be conduct- 
ing only if they intersect the Fermi points. In case that 
{N — M) is divisible by three, there are exactly two such 
conduction bands while otherwise the tube will be semi- 
conducting. Strain and the intrinsic curvature slightly 
shift the Fermi points, which may open a very small 
" semi- metallic" gap, but otherwise leave the electronic 
structure intact i^i^ Due to their large aspect ratio and 
very high structural quality, SWCNTs are often modeled 
as having infinite length and being completely free of de- 
fects. In many situations though, such approximations 
are not entirely justified. 

For example, tubes can be cut into lengths well be- 
low 100 nm allowing for clear experimental observation 
of phenomena related to the one dimensional confine- 
ment of electrons j^iii^i^iiS. Any slight deviations from a 
simple "particle-in-a-box" picture may be due to interac- 
tion effectsiii^ or else must be due to non-trivial bound- 
ary conditions at the ends. Semi-infinite tubes have also 
been shown to have an interesting oscillating electronic 
structure near the ends which even carry some signa- 
tures of interactions J^ii^ However, the tube ends are 
usually not perfect and may contain a number of de- 
fects or can even be closed by structures such as a half 
fullerene, the effect of which we will consider explicitly in 
this paper. Moreover, it is known that even high-quality 
SWCNTs contain on average one defect every 4 jiia}^ 
Imperfections in the bulk such as lattice deformations, 
vacancies or ad-atoms may be induced by strain or irra- 
diation and can drastically modify their electronic and 



transport properties 
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Such irregularities 



can then potentially be used to tailor the properties of 



CNTs2^i2ii2^ for its use in a wide range of applications in 
nano devices In this context, precise understanding 
of the consequences of impurities and structural defects 
becomes especially important. Moreover, the study of 
disorder^^ and single impurities^^i^ in one-dimensional 
systems remains of fundamental interest. It has been 
shown that impurities in general give rise to character- 
istic oscillations of the tunneling density of states that 
reflect the symmetry breaking around the defects^ al- 
beit without considering explicit microscopic models of 
the impurities on the lattice. The amplitudes of density 
oscillations are in turn related to the scattering strength 
and the resulting conductivity.^^ 

In the present work we derive effective lattice models 
of uncoupled chains within the tight binding picture of 
the lowest conduction band. The simple models allow 
the systematic study of the effect from various defects 
in the bulk as well as edge structures on the electronic 
structure of otherwise perfect tubes. In particular, we fo- 
cus on impurity effects in the conduction bands of finite 
metallic SWCNTs with open ends. The non-conducting 
bands of metallic SWCNTs are energetically separated 
by a semiconducting gap that is inversely proportional 
to the radius^ (A ~ on the order of 0.5 eV for 

a radius of 7A. In this work only excitations below this 
energy scale will be considered which are believed to be 
well described by the conduction bands only. The tight 
binding structure of the conduction bands in zig-zag and 
armchair carbon nanotubes can be captured exactly by 
simplified lattice models for SWCNTsi^i^ Our deriva- 
tion results in a lattice of two decoupled chains, which 
allows the incorporation of simple defects or even more 
complicated structures such as caps by adding a small 
number of extra sites. It is demonstrated that, depend- 
ing on the nature of the imperfection under considera- 
tion, the initially independent chains can either remain 
independent or become connected by the presence of the 
defect. Patterns for the local density of states (LDOS) 
near the Fermi level are obtained using such models and 
a straight-forward physical interpretation is possible due 
to the simplicity of the effective theory. At this point it is 
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worth mentioning that such a method accurately repro- 
duces the density patterns for the perfect tube obtained 
by using first-principles calculations in Ref. and the 
results are also consistent with the scanning tunneling 
microscopy (STM) experiments reported in Ref. [9]. As 
such, the simple lattice models here capture the most 
relevant physical effects of defects, but for quantitative 
details future comparison with first principle calculations 
would also be useful. In particular, while our exemplary 
vacancy model only illustrates the first order symme- 
try breaking effect, it would also be interesting to incor- 
porate the relaxation of the structure around defectsi^ 
While this can in principle be considered within the same 
formalism, this extension would require additional input 
from independent ab-initio calculations about the modi- 
fied hopping parameters. 

Electron-electron interaction effects can also be incor- 
porated into the simple lattice models as described in the 
appendix [Al but are ignored for the calculation of the 
STM patterns from defect scattering to first order. This 
is justified because in STM experiments the Coulomb in- 
teraction becomes short range due to screening from the 
substrate and the effective on-site interaction [/gfr = U/N 
is weak for large enough A^j^ It is known that interaction 
effects result in a characteristic amplitude modulation 
of the scattered density wayes^^^^ but appear to 
leave the principle wave pattern intact. Nevertheless, it 
should be emphasized that the models presented here, are 
well suited for the study of strong correlations by apply- 
ing the powerful techniques available for one-dimensional 
systems such as bosonization^ or density matrix renor- 
malization group (DMRG)^^ For the interested reader, 
all additional terms necessary to take electron-electron 
interactions into account in the effective models are listed 
in Appendix [Al 

The paper is organized as follows: In Sec. [Ill we show 
how the effective theory for the perfect CNTs is obtained. 
Then, in Sec. IIIII we incorporate imperfections and caps 
at the tube ends. We proceed similarly in Sec. [IV| but 
now considering structural deformations in the bulk of 
the tube. For each case, we use the effective models to 
obtain electron densities in real space for states near the 
Fermi level. The resulting patterns are the ones expected 
for a typical STM experiment 
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FIG. 1: (Color online) ARMCHAIR: Labels chosen for the 
atoms in the graphene lattice in the case of the armchair CNT. 
The hexagonal network is composed by two triangular sub- 
lattices which are labeled a and p. Within the tight binding 
approximation there is a hopping amplitude t between them. 
The label y is in the direction along the tube axis, and x 
indicates the bonds around the tube. 



suit is a dramatically simplified description of the tubes 
in terms of two chain models. To be able to take advan- 
tage of these simple lattice models, we do not need to 
recur to the widely used linearization around the Fermi 
points j^i^i^i^ The procedure will be outlined for the 
armchair and conducting zig-zag nanotubes for which the 
conduction bands can be easily identified and isolated 
once a Fourier transform around the tube is performed. 



A. Armchair CNT 



II. EFFECTIVE MODELS 

The tight binding band structure of an infinite 
graphene sheet includes a valence band and a conduction 
band that touch at two distinct Fermi points. When this 
graphite monolayer is rolled into a CNT, boundary con- 
ditions are imposed such that only a finite set of bands 
are allowed. As mentioned above, in some cases a subset 
of bands coincides with the Fermi points and the tube 
is metallic. The present section is devoted to show how 
to obtain an effective low energy theory for conducting 
SWCNTs by keeping only its conduction bands. The re- 



Armchair carbon nanotubes are always conducting and 
are defined by a chiral vector with equal indices (A", A").^- 
As usual, the atoms can be divided into two distinct sets 
a and corresponding to the two sublattices of the 
graphene honeycomb structure. It is convenient to la- 
bel each of the a-/3 pairs with a set of coordinates (x, y) 
as shown in Fig. [H The coordinate y is along the tube 
and X is around the tube. The tight binding Hamiltonian 
can now be written as follows 

x=l \y odd 
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FIG. 2: (a) Effective chain lattice model for the perfect arm- 
chair carbon nanotube. The upper chain corresponds to the 
symmetric (S^ mode and the lower one to the ant i- symmetric 
(A) mode. In this case /j. = t. (b) Band structure for the 
conducting modes (effective lattice spacing agff — a/2). 



y even 

"^(^kyf^x^y -^h.C. (1) 



Here a and /3 are the destruction operators in the cor- 
responding sublattice, L is the length of the tube, and 
is the number of bonds around its perimeter. Per- 
forming a partial Fourier transform in the direction x 
around the tube and considering that the only conduct- 
ing modes are the ones corresponding to /c = (one for 
each sublattice),^ the following approximation holds 
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and analogously for /3. Since the k = mode is the only 
one that will be taken into account, we drop the k index 
in the following. Performing the linear transformation 
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the effective conduction band Hamiltonian can be written 
in terms of two completely independent parts, each one of 
them involving only the symmetric (S) or anti-symmetric 
(A) mode 

H,s = Hs + Ha. (4) 
Inserting Eqs. Q and Q into Eq. (HJ we find 

Hs = -tY, (slS^+, + + h.c) (5) 



ky 

Ha = t^(^4^,+i + ^4^, + /i.c.^ (6) 



k. 



Therefore, by keeping only the degrees of freedom of the 
conduction bands, the effective lattice model consists of 
two decoupled chains as depicted in Fig. [21 

Possible more realistic longer range hopping could also 
be included using this formalism. In particular, it is 
also possible to use the exact Wannier orbitals for the 
basis ax^y and Px^y^ e.g. if they are determined from ab- 
initio methods. This in turn would give additional longer 
range hopping terms in the Hamiltonian ([T|) of the con- 
duction band. Since the symmetry is not changed, the 
transformations in Eqs. ([2]) and (|3]) are still equally valid 
and yield again two decoupled chains albeit with corre- 
sponding longer range hopping within each chain. For 
the hybridization due to curvature approximate ax^y and 
Px^y orbitals and the corresponding change in hopping 
have been predicted for the tight binding model^^ which 
are known not to change the symmetry or the electronic 
structure for armchair tubes i^i^i^ 



B. Zigzag CNT 

For the zigzag tube we employ the labeling of the atoms 
shown in Fig.[3l The chiral indices of the wrapping vector 
of the zigzag tubes is given by (A/", 0) where Na is the 
circumference of the tube. Using the notation from Fig. [3] 
the Hamiltonian can be written 
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Introducing the Fourier transform of the fermionic oper- 
ators in the direction around the tube 
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the Hamiltonian becomes 
H = -t 



k I I 
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Only zigzag tubes with N divisible by three are metallic 
and the conduction bands are given by the two wavevec- 
tors k = ±Q = ±^ around the tube.^ Restricting the 
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FIG. 3: (Color online) ZIGZAG: Labeling scheme used for 
the zigzag CNT. Coordinate j marks the sites around the 
tube. The two different sublattices a and /3, now correspond 
to even and odd values of the coordinate / along the tube axis, 
respectively. 



Fourier sum accordingly and using 1 + e^^^s" = s , the 
Hamiltonian for the conduction bands becomes 
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Further simplification is achieved if one introduces the 
transformation 



where 



c±Q,i = e±'^(')C±,, 
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which results in a simple two band nearest neighbor hop- 
ping model 

^eff = E E + ^•^•) (13) 

/i=± I 



By restricting the degrees of freedom to the conduction 
bands, the creation operators on the lattice of the nano- 
tube are therefore reduced to 
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As mentioned above, even and odd chain sites / now cor- 
respond to the two different sublattices of the tube, while 
the two chains correspond to the two allowed wavevec- 
tors k = ±(5 = ±^ around the tube which are exactly 
degenerate. 

For the treatment of defects in the coming sections, it 
will prove useful to introduce the following transforma- 
tion to a parity symmetric basis 
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where is a phase, which is arbitrary for now. It is 
straight-forward to demonstrate that 

^eff = Hg -\- Ha 

I 

and that the creation operators on the carbon network 
are now given by 
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[^1,1 cos(0(/) + Qj - if) 
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sin(0(/) + Qj - if) 



Thus, ip is nothing more than the phase of the wavefunc- 
tion around the nanotube. Since the channels are degen- 
erate, at this point the choice of Lp may seem arbitrary, 
but it is clear that the introduction of structural defects 
and impurities will in general lift this degeneracy. As will 
be illustrated below, in some cases the right choice of if 
allows for the decoupling of the s^p and channels even 
in the presence of complicated structural defects. Note, 
that even small perturbations e.g. coming from the sub- 
strate will in general break the translational invariance 
around the tube, but not necessarily the reflection sym- 
metry. Therefore, in almost all physical situations the 
symmetric/antisymmetric states with some choice of ip 
will be the preferred basis. 

The effective lattice model for this basis and its band 
structure are shown in Fig. [H Interestingly, the effective 
lattice model for the conducting channels has in fact a 
higher translational symmetry than the original lattice. 
The translation by \/3a along the tube in the original 
lattice corresponds to a translation by four sites in the ef- 
fective chains, so that the true wavevector along the tube 
is obtained only after a four- fold zone folding in Fig. [3Kb). 

Rehybridization due to higher bands and curvature 
can also be taken into account within this model for the 
metallic zigzag tubes. The corresponding hopping terms 
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FIG. 4: (a) Effective model for the zigzag CNT. The hopping 
amphtudes are identical for both modes and at half filling 
/X = 0. (b) Since the effective model consists of two identical 
chains, the bands of the s and a channels lie on top of each 
other (effective lattice spacing aeff = aA/3/4). 



have been predicted within the tight binding picture^i^i^ 
and result in a lower symmetry of the model. In partic- 
ular, the effective model in Eq. (p!6|) will also contain a 
very small alternating modulation of the hopping along 
each chain with a corresponding opening of a gap of a 
few meV^^i^i^ We can neglect this effect for excitations 
outside this energy range, which are in fact all but one 
or two for finite tubes of L ^ lOOnm. The degeneracy of 
the two channels is however not lifted by the curvature. 



C. Scanning tunneling microscopy 

The main goal of this study is to predict the effect on 
the electronic structure which can be very well measured 
by scanning tunneling spectroscopy (STS). The corre- 
sponding STM images can be calculated from the simple 
models above in the following way: First the eigenstates 
of finite chains are determined depending on the lengths 
of the tubes, the boundary conditions, and possible de- 
fects, the description of which will be discussed in de- 
tail below. Secondly the resulting eigenstates are trans- 
formed back into the original basis of nanotube orbitals 
{(^x,y^Px,y Cj^i) as described in Eqs. (p1fT8|). Finally, the 
STM images for the eigenstates are calculated by stan- 
dard assumptions of an s-wave orbital on the tip and tt- 
orbitals for the nanotube. For the images presented here 
we use the parameters given in Ref. [4l| in arbitrary units 
for the tunneling current. In order to show the resulting 
impurity effects with a maximum resolution we choose a 
very small distance ^ A of the STM tip to the sample. 
Larger distances naturally result in plots with less res- 
olution but look qualitatively the same. Orbital mixing 
and geometric distortion due to the nanotube curvature 
can also be taken into account ^i^i^ but here the STM 
images are projected onto a flat surface in order to ex- 
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FIG. 5: (Color online) The expected STM image for an arm- 
chair tube with perfectly cut edges. A section near the right 
edge of a tube of total length 30nm is shown, (a) Image cor- 
responding to the excitation of the symmetric mode {S) and 
(b) of the anti-symmetric channel [A), (c) Color-scale of all 
STM images in arbitrary units from to unity. 



plicitly demonstrate the main effects of the defects. The 
procedure above basically determines the LDOS of the 
nanotube by using the square of the eigenstates in the 
effective models. 

In Fig. [5] we show an example for a finite armchair 
tube of 30 nm with sharp edges which corresponds to 
240 lattice sites in each chain of the effective model. The 
symmetric mode for ky = {E ^ 74meV) and the 

antisymmetric mode for ky = {E ^ 93meV) are 

shown. Here the energies are given relative to half- filling 
and we used hopping parameters which correspond to the 
accepted Fermi velocity oi v ^ 8.1 x 10^ The energy 
spacing for a finite tube of 30 nm is about 55 meV for each 
mode (chain). For the armchair tube the symmetric and 
antisymmetric modes are not degenerate. For the zigzag 
tube the two modes are degenerate, but impurities and 
defects generally lift the degeneracy. For example, in a 
30nm zigzag tube the cap described below in Sec. IIIIBI 
lifts the degeneracies by /S.E ~ lOmeV whereas a single 
vacancy on the same tube will have a much smaller effect 
giving /S.E ^ O.lmeV. The color scale given in Fig.[5]has 
been used in all following STM images as well. 

In the STM experiment the LDOS can be measured via 
the differential conductance between the STM tip and the 
nanotube 



\Ei-eV\<5 



(18) 



where 5 is the energy resolution of the measurement. 
For short tubes individual states can be resolved, as was 
demonstrated in several experiments. ^^^^ For longer tubes 
it is still useful to consider individual eigenstates theoret- 
ically, since the resulting experimental signal can always 
be constructed by a simple superposition of the nearly 
degenerate states. 

The impurity effect on the detailed discrete energy 
spectrum is not very easy to measure, while the actual 
modification of the standing waves is much more visible. 
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FIG. 6: (Color online) ARMCHAIR: Expected density pat- 
terns from an STS experiment. The axis of the tube is the 
horizontal direction, (a) Image corresponding to the excita- 
tion of the symmetric mode (S) close to the edge of the nan- 
otube. (b) Excitation of (a) is now phase shifted by = 7r/2. 
(c-d) Same as (a-b) but now for the anti-symmetric channel 
(A). 



which we will discuss in this work. For the pure case 
as shown in Fig. [5] two main features can be observed. 
The short range feature on the atomic scale is charac- 
teristic for the detailed boundary condition as well as 
for the phase of the wavefunction. A long range feature 
shows a rotation of an overall phase of the wavefunc- 
tion, which slowly changes the STM pattern along the 
chain. This long-range feature is not related to impurity 
effects and only depends on the energy relative to the 
half- filled Fermi points. In particular, the rate of change 
with position is always given hy ^ {k — kF)Ay, where 
k-kr ^{E- Ef)/v. 

The short range STM pattern on the other hand 
is characteristic for each mode (symmetric or anti- 
symmetric or a mixture). It is also an indication for the 
detailed boundary condition close to the defects. Typi- 
cal density patterns for the symmetric and antisymmetric 
states are shown in Fig. [6] for an armchair tube (enlarge- 
ment from Fig. [5]). The symmetric states show signifi- 
cant spectral weight on the circumferential bonds, while 
the antisymmetric states always have nodes there. The 
short range patterns are useful to classify the symmetry 
of the states and the overall phase, which only changes 
slowly along the chain as described above. Accordingly, 
the analysis of STM patterns close to boundaries and de- 
fects can be analyzed to classify the nature and symmetry 
of the perturbation as we will see later. 

For the zigzag tube the calculation of the STM images 
proceeds analogously. However, in this case the symmet- 
ric and antisymmetric states for the pure case are degen- 
erate. For most impurities this degeneracy is lifted while 
still preserving the symmetry. However, the choice of the 
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FIG. 7: (Color onhne) ZIGZAG: Density patterns for the 
degenerate s and a channels near Ef- The phase is chosen 
to be (/? = and the tube axis is aligned in the horizontal 
direction, (a) Image corresponding to the excitation of the 
so-mode close to a perfect end of the tube, (b) The state in 
(a) is now dephased by — 7r/2 along the tube, (c-d) Plots 
analogous to (a-b) but now for the ao-channel. 

phase ip in Eqs. ([15]) and (p!8|) is no longer arbitrary and 
depends on the detailed impurity model. Therefore, the 
analysis of the short range STM patterns for the zigzag 
tube gives important additional information. The char- 
acteristic short range patterns of a finite tube are shown 
in Fig. [71 for s^^^o and a^^^o- Other possible values of (f 
will be discussed in the context of defects. 



III. CNT BOUNDARIES 

In this section we examine the consequences of struc- 
tures at the ends of SWCNTs on their electronic prop- 
erties. The tubes can generally have different kinds of 
non-perfect edges such as extra atoms at open ends, pen- 
tagons, or half fullerene caps that can crucially modify 
their behavior. In what follows we study this problem 
in the framework of the chain lattice models derived in 
Sec. [m This will illustrate how simple imperfect ends 
or even complex carbon structures that close the tube 
(caps), can be easily incorporated to the effective Hamil- 
tonians by adding only one or a few sites to the initial 
two chain models representing the clean nanotube. 



A. Armchair 

The simplest non-trivial edge that can be introduced 
is the addition of a single carbon atom at an otherwise 
perfect open end (clean cut) of a SWCNT as shown in 
Fig.[8](a). The extra site can be included in the full hop- 
ping Hamiltonian of an armchair tube in Eq. ([1]) by an 
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FIG. 8: (Color online) ARMCHAIR: (a) Added carbon atom 
at the end of the tube, (b) Extra dimer at the edge. 

additional term 

He = -t + «L2,y=L) c + h.c, (19) 

where c represents the site of the additional carbon atom. 
Since the conduction bands only contain the /c = 
modes, using the approximation in Eq. ([2]) the additional 
term simplifies to 

He.. = {((^Ul + ^Ul)c + h.c.) (20) 

= ^ ( sic ^ h.c) . 

a/TV V ^ J 

Therefore, the electron on the extra site c hops to the 
symmetric /S-mode chain, while the antisymmetric mode 
is unmodified by the additional carbon atom as would be 
expected from the geometry of the problem. The effective 
hopping to the impurity weakens as the radius of the tube 
increases {t' = t/\/N) since the 6'-mode is distributed 
homogeneously around the tube of which only two bonds 
couple to c. 

Another possible defect is the presence of an extra 
bond as shown in Fig. [8Kb). In a similar way as above for 
the additional atom, the extra bond (B) can be included 
into the effective model by adding 

Hb^,, = --^{/3l^Lag+al^j^pB+h.c.) (21) 

— i(a^/3g + h.c.) 

= -J=(^{SiSs-A[A^) + h.c.) (22) 

— t{sis^ — A^^A^). 

As a result the defect is incorporated by adding an extra 
site to both S and A modes with a weakened hopping 
amplitude f = t/\^. It can be easily checked that 



by systematically adding all possible bonds around the 
tube {i.e. extending it), the resulting modification to the 
effective Hamiltonian will be a mere added site, as it 
should be. 

It is worth mentioning that none of the two defects 
considered so far connects the pair of chains of the 
effective model. An example of such mixing produced by 
a defect will be provided later in Sec. IIVB 11 Therefore, 
the modification of the STS spectrum by the extra sites 
can be determined in a straight-forward way: (i) The 
extra carbon atom will have a finite spectral density only 
when a state corresponding to S is excited; otherwise 
it will be completely empty, (ii) STS measurements 
performed over an extra bond are expected to yield the 
same density as the one measured over the bonds at 
the edge of a perfect nanotube that is extended by a/2. 
Otherwise the pattern follows the pure case in Fig. [6] 
with only a slight shift of the phase 0. 

SWCNTs may also be closed at their ends by some car- 
bon structure such as a half fullerene or a nanocone. The 
closure is made possible by the introduction of topologi- 
cal defects (typically pentagons) which induce a bending 
in the carbon network. For the armchair CNT we focus 
here on the (6, 6) nanotube closed by a half fullerene. 
The corresponding cap introduces six pentagon defects 
in order to close the tube, as it can be seen in Fig. [9ja). 

To add the caps to the effective model, the only nec- 
essary modification is the inclusion of two extra sites at 
each end of the -S-mode chain. The first one of them rep- 
resents the k = mode of the ring formed by the atoms 
at the tip of the 6 pentagon defects near the end as shown 
in Fig. [91(a). The last site corresponds to the k = mode 
of the hexagon at the top of the cap. The values of the 
hopping and chemical potential for each one of these ex- 
tra sites will be different from the bulk values t' = \/2t^ 
fi' = 0, t^^ = t, fi" = 2t as indicated in Fig. [9l Due to 
the pentagon defects of the cap, the chain corresponding 
to the antisymmetric mode is not modified by the addi- 
tion of this structures at the ends of the tube. It is easy 
to understand this by considering that the atoms at the 
tips of the pentagon defects are added in the same way 
as explained in Eq. ([2Q|) for a single pentagon. In this 
picture the A-mode becomes then completely decoupled 
from the half fullerene cap.^^ As a result, the effective 
lattice model that includes the cap looks like the one 
depicted in Fig. [H^b). Although we have here outlined 
the procedure for the specific cap shown in Fig. [9](a), 
the same basic ideas can be applied to construct effective 
models for nanotubes of larger radius. Capping struc- 
tures for such tubes will in general be composed by a 
larger number of carbon atoms; therefore, the resulting 
effective models will require more additional sites for its 
description. 

For the individual states near a capped nanotube very 
clear effects in an STS experiment are expected. Namely, 
if the applied voltage matches the energy of one of the 
states corresponding to the A-mode, the cap will appear 
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FIG. 9: (Color online) ARMCHAIR: (a) Half fullerene cap at 
the end of a (6, 6) carbon nanotube. (b) Low energy effective 
model for the (6, 6) nanotube including a half fullerene cap. 
The sites corresponding to the cap are shown in red (gray). 
The A- mode is not modified by the cap. Parameters for the 
added sites are: t' = V2t, ^' = 0, t" = t, = 2t. (c) Elec- 
tronic density for an eigenstate of the symmetric {S) mode, 
obtained using the effective model in (b). For illustration pur- 
poses sites corresponding to the cap have been flattened out. 
(d) Same as in (c) but now for the A-mode; as expected, the 
cap is empty in this case. 



completely empty. On the other hand, exciting a sym- 
metric state would be clearly recognizable by the appear- 
ance of a finite amount of electron density distributed 
over the cap. Our density maps (simulated STM images) 
clearly illustrate this effect as shown in Fig.[9](c) and (d). 
The density over the body of the tube is distributed very 
similar to the pure case in Fig. [6l 





FIG. 10: (Color onhne) ZIGZAG: (a) A vacancy at the edge 
of the tube, (b) Extra dimer producing a pentagon defect at 
one end of the CNT. 



their corresponding low energy models are as simple as 
possible. 

We first consider the absence of a carbon atom at the 
end as illustrated in Fig. fTOlf a). This sort of imperfection 
involves the elimination of a pair of bonds that joined the 
missing atom at site / = 1 with a pair of atoms at / = 2 
plus the addition of an extra hopping matrix between 
them. The resulting effective impurity Hamiltonian for 
such a vacancy (V) is 



2t 

N 



Sal 



(23) 



where we have chosen cp = 7r/6. The first term corre- 
sponds to the elimination of the atom and the second to 
the inclusion of the extra bond shown in Fig. [TOlfa). 

Another possible defect is the addition of an extra bond 
as the one shown in Fig. [TOlfb). In the effective models 
it is represented by couplings of both chains to the bond 
sites, corresponding to the impurity Hamiltonian 



Hb,,, = t (blb_ - 6^6+) 



^0,1^+ 



I.e.) , (24) 



where, 



B. Zigzag 

As already mentioned in Sec. [Ill B there is an arbitrary 
phase cp that defines the transformation to the effective 
two chain models. In what follows we will choose cp for 
each structural defect that we treat in such a way that 



bl 



b\±bl 
V2 ■ 



(25) 



and bi^2 are the sites of the extra bond. We note that the 
choice cp = keeps the two modes completely decoupled 
even in the presence of a bond defect at the edge. There- 
fore, the STS signal will look very similar as in the pure 
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FIG. 11: (Color online) ZIGZAG: (a) Half fullerene cap at 
the end of a (9, 0) CNT. (b) Low energy effective model for 
the conducting zigzag nanotube including the half fullerene 
cap with parameters as defined in Eq. (|27)) . The last site of 
the s (a) chain corresponds to the /c = (A; = tt) mode of the 
hexagon located at the top of the cap structure, (c) Local 
density of states for an sil state and (d) for an aiL state. 



case in Fig. [7] (except for a small change of the phase 
along the chain). 

To demonstrate the inclusion of a structure that closes 
the tube, we focus on the case of the (9, 0) nanotube de- 
picted in Fig. [TlT a). Such a half fullerene cap can be 
constructed starting from the perfect CNT in the follow- 
ing way: First, three vacancies are introduced at one end 
(/ = 1), each one of them producing a pentagon defect as 
the one discussed above for the vacancy. Then, each one 
of the 6 remaining atoms at the end is linked to the atoms 
that form the hexagon at the top of the cap. Choosing 
again cp = ^ it turns out that each one of the afore- 
mentioned "cap vacancies" (CV) correspond exactly to a 
term like the one obtained above, in such a way that the 
part of the cap Hamiltonian corresponding to all three of 
them is 



(26) 



Furthermore, it can be shown that the S2l and ajL chan- 

' 6 6 

nels of the effective model couple only to the k = and 
k = TT modes {dl ^) of the hexagon at the top of the 
cap. The resulting effective Hamiltonian for the cap of 
the (9, 0) nanotube becomes 



HcAP = Hc\ 



where, 



2i(44 - 4do) (27) 



dk 



1 



(28) 



and dj are the sites of the hexagon at the top of the cap. 
By conveniently choosing Lp = 7r/6 we have again man- 
aged to express the Hamiltonian in terms of two com- 
pletely decoupled modes even when the cap is present. 
The resulting model requires the addition of a couple of 
sites and a local change in the values of the parameters 
as it is illustrated in Fig. [TTlfb). As for the armchair 
nanotube, we have described here the procedure only for 
one specific cap. Nevertheless, the same scheme can be 
followed to model caps of tubes of larger radius that will 
naturally require a greater number of extra sites in the 
effective model. 

The LDOS around the cap of the (9,0) nanotube is 
displayed in Fig. \TT\ c) and (d) for two different energies 
near the Fermi level. In this case both s and a states have 
a finite density at the cap. Notice that due to the cp = 7t/6 
phase shift the forms of the patterns are inverted relative 
to Fig. [71 the ones of S2l (aiL) look like the ones of ao 
[So). This can be easily understood from Eq. ([T5|) since 
the shift ofip — ir/Q together with a translation of one site 
around the tube converts the cosine and sine functions 
into one another. The effect is clear when comparing 
Fig. with Fig.Hl 



IV. DEFECTS IN THE BULK 

The honeycomb structure composing the CNT may 
present a variety of structural defects which are known 
to significantly modify its electronic properties J^i^ It is 
possible to model any kind of defect in a straight-forward 
way using the chain lattice models introduced in Sec. [TTl 
As an example we will show explicitly how to include two 
prototypical impurity models: the Stone- Wales and the 
single vacancy defects as shown in Fig. [121 We first show 
how the low energy models look once the defects are in- 
cluded, and then use them to obtain electron densities for 
the states closest to the Fermi energy. As already men- 
tioned in the introduction, relaxation^ will in general 



10 




(a) armchair (b) 




armchair 



FIG. 12: (Color online) Two types of defects in the bulk of car- 
bon nanotubes. The arrows indicate the direction of the axis 
for the two classes of tubes, (a) Stone- Wales defect formed 
by a pair of pent agon- heptagon deformations. Defective sites 
are depicted in red and they correspond exactly to the ones 
added to the effective model, (b) Vacancy defect produced 
by taking away a single carbon atom, breaking the a-p sub- 
lattice symmetry. The big dot (blue) represents the missing 
atom also in the following figures. 



also occur around defects which cannot be implemented 
in the formalism without additional input from ab-initio 
calculations. Nonetheless, the examples we consider here 
demonstrate the general method also for more compli- 
cated cases, and also show the most relevant effects that 
are expected from impurities of this symmetry class. 



A. Stone- Wales defect 

One of the simplest deformations that can be in- 
troduced to the perfect graphene lattice is a pair 
of pentagon-heptagon defects as the one shown in 
Fig. [Ha). These are usually called Stone- Wales"^"^ (SW) 
or 5-7-7-5 defects and can be thought of as a simple 90° 
rotation of a single bond. Stone- Wales defects have low 
energy and they do not modify the helicity of the tube. 
They are also known to be responsible for the initiation 
of the plastic deformation of CNTs as the bond aligns in 
the direction of the applied strain4^ 



1. Armchair 

Incorporation of the 5-7-7-5 defect to the effective 
model of the armchair nanotube is achieved by adding 
a couple of sites that are connected to the 5'-mode chain 
with the parameters = t{N - l)/N, ii' = fi{N - 
1)/N, t" = as shown in Fig. [IHa). This can 

be derived in two steps: First, eliminating one bond has 
the effect of locally weakening the hopping and chemical 
potential of the effective chain sites. Then, the rotated 
bond is included by adding a couple of extra sites corre- 
sponding to the atoms of the defective bond. As pointed 
out above, the pentagon defects along the tube are not 
seen by the antisymmetric mode; therefore, the second 
step does not modify the A chain. 

Again, there are interesting effects that can be ob- 
served in simulated STM images. Tuning the voltage to 
an energy corresponding to an eigenstate of the A chan- 
nel, these sites will be significantly occupied as seen in 
Fig. [T3lf b). Otherwise, if the applied voltage matches the 
energy of the antisymmetric mode, the defective bond 
will appear as completely empty as shown in Fig. [T3lfc). 
In both cases, away from the defect the density pattern 
will look essentially identical to the one of the perfect 
tube as shown in Fig. [6l 

2. Zigzag 

The situation for the zigzag tube is rather different. 
First we choose = 7r/6 for the modes and in the 
transformation ([15]), since it yields the simplest possible 
low energy model. In such a basis the elimination of a 
bond translates into a weakening of a few hopping ele- 
ments of the s chain but leaves a intact. Furthermore, 
the subsequent addition of the defective bond involves 
connecting the s (a) channel to the 6+ (6_) mode of the 
added bond, where h± is defined by Eq. ([25]) . The re- 
sulting form of the effective model with the parameters 
t' = t{N - 1)/7V, t+ = t/\^ and t- = ±iA/3t+ is shown 
in Fig.[T4lfa). 

Patterns obtained from an STM experiment around 
the imperfection for energies close to the Fermi level are 
expected to look like Fig. [TWb) and (c). Depending on 
the symmetry of the state and the position of the defect, 
the Stone- Wales bond can be observed to accumulate a 
large amount of density as in Fig. [T4lfb) or in other cases 
become completely empty as in Fig.[lU^c). Unlike for the 
armchair case, in the zigzag CNT the lack of density on 
the defective bond is not an exclusive property of a state 
with a certain symmetry. 

B. Vacancy defect 

Irradiation of a carbon nanotube by Ar+ ions may re- 
sult in the presence of metastable but long-lived single va- 
cancy defects on its surface^ as illustrated in Fig. [T2lf b). 




FIG. 13: (Color online) ARMCHAIR: (a) Effective lattice 
model including the Stone- Wales defect. The extra sites (in 
red) correspond to the two carbon atoms of the defective 
bond. Some hopping amplitudes and chemical potentials are 
modified around the deformation: t' = t{N — l)/N, ii' = 
fi{N - 1)/N, t" = t^/2/N. (b) Electron density around the 
defect for the S state at the Fermi level. As expected, there is 
a finite occupation of the defective bond, (c) Same as in (b) 
but now for the A state. The defective bond is now completely 
empty. 



The absence of an atom in the carbon network obviously 
does not change the overall chiral vector of the tube, but 
it will still have some pronounced consequences on its 
electronic properties. For example, recent theoretical^i 
and experimenta l^^i^^ works have shown that nanotubes 
enter the strong Anderson localization regime when the 
concentration of defects is large enough. 



1. Armchair 

To include the vacancy in the effective chain model a 
single extra site must be added, which corresponds to the 
remaining atom of the defective bond with coupling pa- 
rameters = t{N-l)/N, ii' = /i(7V-l)/iV, t" = t/^/2N 
as shown in Fig. [T5lfa). It is interesting to note that this 
time the S and A chains become connected as a con- 
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FIG. 14: (Color online) ZIGZAG: (a) Effective lattice model 
for the zigzag tube including the Stone- Wales defect. Addi- 
tional sites correspond to the symmetric (6+) and antisym- 
metric {h-) modes of the defective bond. The values for the 
locally modified parameters are t' = t{N — l)/N, t+ = t/VN 
and t- = ±i>/3t+. The hopping direction in which t_ has a 
positive sign is indicated by arrows to b-. (b) Electron den- 
sity around the defect for a sil state close to the Fermi level, 
(c) Same as in (b) but now for a aiL state. 



sequence of the introduction of the defect. Therefore, 
eigenstates will not necessarily correspond to the elec- 
tron occupying a single mode, but they will in general be 
a mixture of both. It turns out that there are now three 
types of possible states: S'-like, A-like and mixed (M). 
The first two correspond to states in which the electron 
occupies primarily only one of the original channels as 
illustrated in Fig. [T5lfb) and (c). Mixed states are also 
possible, where the eigenstate is a linear combination of 
the S and A modes. It is found that the M states com- 
pletely break the sub-lattice symmetry of the hexagonal 
lattice, resulting in an accumulation of the electron den- 
sity on the sublattice opposite to the one of the vacancy 
as shown in Fig. fTSlf d). It appears that the effective extra 
site enables a resonance at certain energies. A related ef- 
fect of a strongly enhanced spectral weight around vacan- 
cies at certain energies has been predicted for graphene 
using a linearized band structure4^ 

For quantitative estimates more realistic microscopic 
models must also take relaxation effects into account 
which are especially important for vacancy defects. This 
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FIG. 15: (Color online) ARMCHAIR: (a) The vacancy is 
included into the low energy effective lattice model. An ex- 
tra site represents the nearest neighbor of the vacancy in the 
direction perpendicular to the tube's axis. Parameters are 
modified around the imperfection: t' — t(N — 1)/N, /j,' — 
fi{N - 1)/N, t" = t/V^N. (b) Electron density around the 
defect for an eigenstate which looks basically exactly the same 
as one corresponding to the pure symmetric channel (iS-like). 
(c) Same as in (b) but now the density looks as the one of 
the A mode (A- like), (d) There is a third group of eigenstates 
that combine both channels in such a way that the density is 
accumulated in one of the sublattices; this is a manifestation 
of the broken sublattice symmetry by the vacancy defect. 




FIG. 16: (Color online) ZIGZAG: (a) Chain lattice model 
including the vacancy defect. There is a slight weakening of 
the hopping along the S|. channel from t to t' — t{N — l)/N. 
(b) Local density of states around the imperfection for a level 
corresponding to the chain, (c) Density pattern around 
the vacancy for an aiL state. 



will generally modify the couplings and STM images in 
a finite range, but the resulting eigenstates can be classi- 
fied into one of the three types (A), (S), or (M) with the 
characteristic general properties as shown in Fig. [T5l 



2. Zig-zag 

For the case of the zigzag nanotube, the vacancy does 
not require an additional site in the effective model since 
it represents only a slight modification on one effective 
chain lattice site. In this case, the choice (/? = 7r/6 again 
yields the simplest effective model. In fact only the sz. 
mode is modified by the presence of a vacancy which 
causes a slight weakening of local hopping terms t' = 
t{N—l)/N as shown in Fig. fTHlf a). We note that the same 
impurity terms were also present when we considered the 
Stone- Wales defect, which is not surprising given that it 
is nothing more than two consecutive vacancies plus the 
attachment of the defective bond. 

As seen in Fig. \Wl h) and (c), the broken sublattice 
symmetry is not manifest in the density patterns. Nev- 
ertheless, the presence of the vacancy has the effect of 
pinning the phase of the wavefunction around the nan- 
otube. Remarkably, the antisymmetric channel is not at 
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all affected by the vacancy, except of course for the miss- 
ing spectral weight at the vacancy. As discussed before, 
the choice of = 7r/6 makes the antisymmetric state ap- 
pear like the cp = symmetric state in Fig. [71 and vice 
versa. 

As mentioned above relaxation is important also in the 
zig-zag case. However, since the parity symmetry of the 
system is not broken by a vacancy in the zig-zag tube the 
channels are bound to remain decoupled and the general 
classification in Fig. [3 still holds. 



V. SUMMARY AND OUTLOOK 

In summary, we have shown how defects and capping 
structures can be added to effective two chain models 
of conducting SWCNTs. We have demonstrated the 
method by incorporating half fullerene caps and Stone- 
Wales and vacancy defects to both the armchair and 
conducting zigzag nanotubes. Using the resulting chain 
lattice models of the conduction bands we obtained the 
LDOS around the imperfections for states close to the 
Fermi level. 

Besides their simplicity the most remarkable feature of 
the effective models is the fact that the two chains are 
completely decoupled and also remain decoupled even in 
the presence of most defects and additional structures at 
the end considered here. Especially in the zigzag case, 
the defects in fact dictate the proper choice of the ba- 
sis, which then represent the two independent channels. 
Surprisingly, a single vacancy leaves all states in the an- 
tisymmetric channel of metallic zigzag CNTs completely 
unchanged. Therefore, transport in that channel should 
be possible without scattering, which may have strong 
consequences for the resulting conductivity. However, 
this reduction of scattering is only possible if vacancies 
are spaced so that they do not affect each other and if 
there are no other strong perturbations which would pin 
the phase (p around the tube. Interactions will introduce 
a coupling between the independent chains as discussed 
in the appendix, but of course even in that case it is 
useful to start in the basis of decoupled chains. 

At this point it should be emphasized that the meth- 
ods presented above are not restricted to the examples 
explicitly shown here, but the transformations in Sec. [Ill 
give a quite general framework which allows for the treat- 
ment of a wide range of possible imperfections. For ex- 
ample, the procedure illustrated in the present paper can 
be extended to tackle problems such as junctions of (con- 
ducting) CNTs of different chirality, presence of magnetic 
impurities, or transport related problems, just to mention 
a few interesting prospects. 

From the point of view of applications, it is clear that 
this approach can be useful to assess the effect that a 
particular deformation or impurity has on the electronic 
properties of nanotubes in a very simple and precise man- 
ner. It is thus a promising tool to be employed for the 
design of CNT based microelectronic devices. 



On the other hand, from a more fundamental perspec- 
tive, since the resulting models are one-dimensional, it 
will be certainly interesting to take into account strong 
correlation effects in future works by using the techniques 
available for such systems {i.e. bosonization or DMRG). 
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APPENDIX A: INTERACTIONS 

In the main text we have neglected interactions be- 
tween electrons and worked exclusively with the result- 
ing free electron models. We will now show that in- 
cluding interactions into the effective Hamiltonian can 
be done in a straight-forward way. Although in princi- 
ple we could consider any electron-electron potential, for 
definiteness we will treat only on-site density-density in- 
teractions (U). To illustrate how a defect modifies the 
interacting term, we will explicitly show the case of the 
Stone- Wales deformation. 



1. Armchair 

For the armchair carbon nanotube an on-site interac- 
tion is given by 

= UJ2 (4,t"-- T^i"--! + 4t/3.- 1/3^4/3?, i) (Al) 

r 

Using Eq. ([2]) the conduction band part becomes, 
y 

In terms of the S and A modes this term looks slightly 
more complicated, 

y 

+ ^sAivWsAiy) ^^Asiy)^As(y) 
+ ^ssiy)^AAiy) ^^AAivWssiv) 

+ ^SAiV^Asiv) + ^AsivhsAiy)) ^ 

where, 

n^PQ = PlQ. (A4) 

As we can see, even though rewriting the Hamiltonian in 
terms of the S and A modes simplifies its non interacting 
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part, the interaction terms become more convoluted. It 
is interesting to note that the effective on-site interaction 
becomes inversely proportional to the number of bonds 
around the tube (/7eff = U/AN). Such a dependence is 
not surprising: the occupation of a site in one of the 
effective low energy channels represents the presence of 
a single electron distributed around the tube's N bonds. 

To include the Stone- Wales defect we first note that a 
(normal) bond around the nanotube is missing and there 
is thus a slight decrease of the on-site interaction at the 
corresponding site, 

= ^^^^C/eff. (A5) 

Furthermore, we have to add the interaction terms that 
act on the added (defective) bond, 

Hr = U{nlnl+nlnl). (A6) 

Remarkably, while the interaction in the bulk of the tube 
is reduced by a factor of 1/N in Eq. (jA3p . the effective 
interaction on the Stone- Wales bond is not rescaled and 
therefore is bound to play a much stronger role. 



where, 

= Cl,C,,a. (A9) 

We find that the effective interaction is again inversely 
proportional to the tube's circumference (/7eff = U/N). 
Note also that the initially independent chains of (p!3|) 
become correlated by the scattering terms that involve 
both + and — modes. 

The first effect that the Stone- Wales defect will have 
on this model is a weakening of the interaction strength 
(t/eff) similar to the one derived for the armchair tube, 
but now it affects two contiguous sites, and + 1 (see 
Fig.[l4j). They wih both have. 



UeE = ^^^t/eff. (AlO) 

Furthermore, a couple of terms of order 0{N~'^) arise, 
which contain local interaction vertices of a kind that 
was not present for the pure system. 



2. Zigzag 

The interaction Hamiltonian for the zigzag CNT is, 

^'"' = ^ESvT^^VT4,U^^Vi- (AT) 

Using the two band approximation and transformation 
in Eq. ([T4j) the effective interaction term becomes. 



^eff =^E (nU(Oni+(0+nI_(Oni_(0 (A8) 



+ n^+(0ni_(0 + nI_(0n^+(0 
+ n;_(0ni+(0 + nl+(0ni_(0), 
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one for each of the defective sites mentioned before (/ = 

The last and most trivial part of the interaction term 
coming from the impurity corresponds to the on-site 
terms of the defective bond. Such a term will be ex- 
actly the same as the one given for the armchair carbon 
nanotube in ()A6p . again with a not rescaled stronger in- 
teraction strength. 
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